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1. Introduction 

The polarizabilities of a composite system such as the pion are elementary structure constants, 
just as its size and shape. They can be studied by applying electromagnetic fields to the system, that 
is, by the Compton scattering process 7+ % — > 7+ % or the crossed-channel reactions 7+ 7 <^ % + n. 
The low-energy theorem expresses the zeroth- and first-order terms of the amplitude by the charge 
and the mass of the pion. The second-order terms in the photon energy describe the response of 
the pion's internal structure to an external electric or magnetic dipole field, they are proportional to 
the electric (a) and magnetic (j8) dipole polarizabilities. In this contribution we concentrate on the 
polarizabilities determined by forward (a + j8) and backward (a — j8) scattering. 
Within the framework of the partially conserved axial-vector (PCAC) hypothesis and current alge- 
bra, the polarizabilities of the charged pion were related to the radiative decay n + — ► e + v e y [jl|]. 
Chiral perturbation theory (ChPT) at leading non-trivial order, 0(p 4 ), confirmed this result, a K + = 
—fin+ ~ @L where l& = (4 — fs) is a linear combination of scale-independent parameters of 
the Gasser and Leutwyler Lagrangian [gj. At &{p A ) this combination is related to the ratio of 
the pion axial-vector form factor F A and the vector form factor F v of radiative pion beta decay, 
Fa/Fv = 1a/ 6 [||]. Once this ratio is known, chiral symmetry makes an absolute prediction at 
0(p A ), tt n + = 2.64 ±0.09, here and in the following in units of 10~ 4 fm 3 . Corrections to this 
leading-order result were calculated at ^(/? 6 ) and turned out to be rather small [Q, ||]. This makes 
the following predictions for the polarizabilities a very significant test of ChPT [j|] : 

= 0.16±0.1, (1.1) 
a n+ -fi n+ = 5.7 ±1.0. (1.2) 

The results of ChPT are in sharp contrast with recent predictions based on dispersion relations 
(DRs) [§, 

a n+ = 13.60±2.15. (1.3) 

In this work, the dispersion integrals are saturated by various meson contributions in the s and t 
channels. The free parameters are fixed by the masses, total widths, and partial decay widths of 
these mesons at resonance. However, the extrapolation to energies below and above resonance is 
performed with specific resonance shapes whose analytic properties leave room for considerable 
model dependence. 

In the present contribution we address the conflicting results obtained by ChPT and DRs. Section || 
gives a brief introduction to the Mandelstam plane and the kinematics used to describe the scatter- 
ing amplitudes. In Sec. || we present the elements of the DRs used by Fil'kov and Kashevarov ^ 
and investigate the analytic structure of this model along the lines of our earlier work [R]. Our 
conclusion in Sec. [| summarizes the arguments against the predictions of Ref. 

2. Kinematics and Mandelstam Plane 

The Mandelstam variables for Compton scattering, y{k) + 7t(p) — > y{k') + Ti(p'), 

s=(k + p) 2 , t = (k-k') 2 , u=(k-p') 2 , (2.1) 
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Figure 1: The Mandelstam plane for Compton scattering on the pion. Left panel: The hatched areas for 
s > m 2 , t > 4m 2 , and u> m 2 show the physical regions in the s, t, and u channels, respectively. Right panel: 
The hatched area forms the triangle bounded by the thresholds of two-pion production, s = Am 2 , t = Am 2 , 
and u = Am 2 . See text for further explanation. 

are constrained by s + 1 + u = 2m 2 , where m is the pion mass. The crossing-symmetric variable v 
is defined by 

s — u 

V = — . (2.2) 



Am 



The two Lorentz-invariant variables v and t span the Mandelstam plane shown in Fig. [I]. They are 
related to the initial (Ey) and final (EL) photon lab energies and to the lab scattering angle 6 by 



V=Ey + - = -(E Y+ E'y), 
t = -4EyE' y sin 2 (6/2) = -2m(Ey-E' y ) 



(2-3) 



The scattering matrix can be expressed by 2 independent amplitudes, denoted by M ++ (v,t) and 
M + ~(v,t), the superscripts indicating the helicity of the photons in the channel 7 + 7— > % + %. 
Because of the crossing symmetry they satisfy the relation M +± (— v,t) = M +± (v,t). 

The Mandelstam plane is displayed in Fig. [j] The hatched areas in the left panel show the 
physical regions in the s, t, and u channels. The dashed lines mark the onset of inelasticities due 
to the thresholds of two-pion states with invariant mass > 4m 2 in the respective channels. The 
polarizabilities are obtained in the quasi-static limit, that is, in the limit of approaching the origin 
of the Mandelstam plane, s = u = m 2 or v = t = (solid circle). The polarizabilities may be 
calculated by dispersion integrals running along lines through the origin, e.g., v = 0, t = 0, and 
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u = m 2 . The right panel of Fig. [j] highlights the triangle below the thresholds: within the dashed 
area, the non-Born scattering amplitudes for real t and v are real functions without any singularities. 
Therefore, these amplitudes can be Taylor expanded about the origin of the Mandelstam plane, 
with a convergence radius given by the border lines of the triangle. The leading terms of these 
Taylor series describe Thomson scattering, the polarizabilities can be read off the subleading terms. 
Contrary to this basic requirement, the model of Ref. ^ contains lines of square-root singularities 
passing through the internal triangle. The vector mesons are modeled with a factor 1 /^s, which 
also leads to a factor 1 / yfu by the crossing symmetry, and the scalar a meson contains a factor l/y/i 
leading to a divergence at the origin of the Mandelstam plane, at which point the polarizabilities are 
determined. In the following we concentrate on the vector mesons, in particular the p(770). The 
case of the a meson has been discussed in Ref. P], most of the other intermediate meson states in 
the work of Ref. [§] have similar singularities. 



3. Model of Fil'kov and Kashevarov for Vector Mesons 

The model of Ref. ^ describes the contribution of a vector meson V = {p,co} to Compton 
scattering by an energy-dependent coupling constant g(s) at the vertex y+K ^ V and a vector 
meson propagator 1/{(M — iF(s)/2) 2 —s}, with M the mass of the vector meson and T(s) its 
energy-dependent width. The term quadratic in F(s) is neglected (small-width approximation). In 
the notation of Ref. [0], the vector meson contributions to the amplitudes take the form 

The width T(s) has the energy dependence of a P wave at threshold, 

4,„ 2 



r(s > = U^wJ r »- <3 - 2) 

with To the width of the vector meson V at resonance, s = M 2 . Furthermore, the square of the 
energy-dependent coupling constant, g(s) 2 , is defined with an s -1 / 2 singularity in order to obtain 
convergence of the dispersion integrals at large energies, 

6kM ( M s? " 



with Ty the partial decay width for V — > % + y. 

The vector meson contributions to the polarizabilities are derived from the amplitudes by 

a + B = ^-M + -(s = m 2 ), a-B = -J— M ++ (s = m 2 ). (3.4) 



Combining these equations with Eq. (|3JJ), we obtain 

a + B = -(a-B) =► a = 0. (3.5) 

The (quasi-static) electromagnetic transition from the pion (J p = 0~) to the intermediate vector 
meson (J p = 1~) is a magnetic dipole transition yielding a paramagnetic contribution to B and, as 
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meson 


a + P a-fi 


a J8 


P 
00 


0.063 -1.15 
0.721 -12.56 


-0.54 0.61 
-5.92 6.64 



Table 1: Results of Ref. []6[] for the contributions of p and a mesons to the polarizabilities of charged and 
neutral pions, respectively, in units of 10~ 4 fm 3 . 
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Figure 2: Dispersion relations for Compton scattering on the pion. Left panel: The s channel singularities 
introduced by Ref. [^] are the vector meson resonance near s = M 2 , the onset of the physical cut at s = 4m 2 , 
and the onset of the unphysical cut at s = 0. In order to determine the polarizability at (v = 0,f = 0) (solid 
circle), the dispersion integral can be obtained along different paths, e.g., t = 0, u = m 2 , and v = 0. Right 
panel: Contour integrals in the complex v plane, related to forward DR (path along t = 0). On provision that 
the physical Riemann sheet is free of singularities except for the physical cuts describing particle production, 
the Cauchy integral along the small circle about v = can be replaced by the dispersion integrals over 
the discontinuities along the cuts (dashed lines) and a big circle whose contribution should vanish for an 
infinitely large radius. However, the model of Ref. [^] contains additional square-root singularities at s = 
and u = 0. See text for further explanation. 



a consequence, a ratio R= (a — j3)/(a + /$) = — 1. To the contrary, Fil'kov and Kashevarov predict 
a ratio R as —20 and a large electric polarizability a (see Table [T]). Even more surprising, the latter 
carries a negative sign, a result possible only in a relativistic quantum field theory. In conclusion, 
the results of Ref. ^ are at variance with Eq. (3.5). 

The discussed discrepancy is related to the unphysical singularities at s = 0, t = 0, and u = 0, very 
close to the origin of the Mandelstam plane at which the polarizabilities are determined. The left 
panel of Fig. @ displays the singularities of the model in the s channel, the unphysical singularity at 
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s = 0, the threshold for two-pion production at s = 4m 2 , and the vector-meson pole position near 
s = M 2 . Also shown are the paths of integration for three DRs, along t = 0, u = m 2 , and v = 0. 
Fil'kov and Kashevarov evaluate the dispersion integral over the physical cut, from s = 4m 2 to 
infinity, they ignore the s~ l l 2 singularity and the unphysical cut resulting from it. However, the 
factors s~ l l 2 and m -1 / 2 are not small details but necessary to fulfill the premises of Titchmarsh's 
theorem, the square integrability of the amplitude along any line parallel to the real axis. Moreover, 
even for a model fitting the data in the physical region, one should not trust the extrapolation to the 
origin of the Mandelstam plane if the model has a singularity at or close to this point. Setting the 
imaginary parts of the amplitudes equal to zero for Re[s] < 4m 2 , Fil'kov and Kashevarov tacitly 
annul Titchmarsh's theorem and at the same time create a non-analytic amplitude, for which no 
DRs exist. The right panel of Fig. |2] illustrates the problem in the complex v-plane. Let us start 
from a contour integral along a small circle about v = 0. If the function F(v) is analytic within the 
circle, the function at v = is given by Cauchy's integral, 

Assuming that the physical sheet contains only the singularities of two-pion production at v = 
iVthr = ±| m, we can blow up the small circle to the big contour in the right panel of Fig. ^, with- 
out changing the value of the integral, F(0). On condition that the function is square integrable, the 
contribution of the large circle (solid part of the big contour) vanishes if its radius goes to infinity. 
The function F(0) is then determined by the integrals over the discontinuity of Im[F(v)] on the 
physical cuts — oo < Re(v) < — Vthr and v t h r < Re(v) < °° (dashed part of the big contour). Because 
the discussed model has additional square-root singularities at v = ±Am, additional unphysical 
cuts appear, which we suggest to draw along the real axis such that at least the neighborhood of 
V = remains an analytic region. With these provisions Titchmarsh's theorem assures that DRs 
exist, however at the expense of introducing unphysical cuts. 

Although the square-root factors are essential to provide the required convergence at infinity, 
Fil'kov and Kashevarov neglect the large contributions from the associate discontinuities. As an 
example, such neglect could be turned into the statement Im[F(v)] = if — v t h r < Re[v] < v t hr- 
The result would be two "walls" in the complex plane with a non-analytic strip in between. Quite 
generally, any statement that the imaginary part of an analytic function be zero in a certain area 
makes this area (and therefore the function itself) a non-analytic one, because an analytic function 
can only contain isolated singular points. 

The analytic structure of the model amplitudes is further illustrated in Fig. [| The left panel shows 
the contour plot of Im[M + ~ (s)] in the complex plane, s = x + iy. The physical (right-hand) cut with 
a maximum discontinuity near x = M 2 « 0.55GeV 2 is clearly seen. However, the unphysical (left- 
hand) cut ranging from s = to the left turns out to be of similar importance. Near x = — 1 1 GeV 2 , 
we also observe a "bound state" embedded on the left cut. This unphysical pole is a consequence 



of the energy-dependent width in the resonance propagator of Eq. ( |3.1[ ). The forward polarizability 
a + /3 is determined at the point (x = m 2 , y = 0), squeezed in between the two cuts. The right 
panel of Fig. || shows the same plot for Im[M ++ (.s)], which is related to the backward polarizability 
a — f5. Because of the additional factor s in this amplitude, the p -meson resonance appears sup- 
pressed. 
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Figure 3: Contour plots for the imaginary parts of the p meson contributions to M + (s) (left panel) and 
M ++ (s) (right panel) in the complex plane with coordinates s = x + iy, all in units of Ge V 2 . 



Integrand for a+/3 Integrand for a-fi 




Figure 4: The integrands of the dispersion integrals for the polarizabilities a + j3 (left panel) and a — /3 
(right panel), as function of x in units of GeV 2 . See text for further explanation. 



Figure ^displays the integrands of the dispersion integrals for the polarizabilities a ± j3 as function 
of x = Re [s]. The integration yields the following ^-channel contributions for the p meson: 

a + p = 0. 17 (I.e.) + 0.03 (r.c.) = +0.20 , 

a-fi = 0.98 (I.e.) -1.18 (r.c.) = -0.20, (3.7) 
where I.e. and r.c. denote the contributions of the left and right cuts, respectively. In agreement 



with Eq. (3.5), the total result yields a paramagnetic contribution of the vector meson to the po- 



larizability. In particular, R = (a — j8)/(a + j8) = —1, whereas the right-hand contributions lead 



to a ratio close to R = —40. We note that the values for a + j8 given in Eq. (3.7) include only the 
s-channel contribution of the p meson, whereas the predictions of Ref. [Q| listed in Table [I] contain 
both s and u channels. 
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a + /3 




a — 


P 




model 


real 


r.c. 


l.c. 


rest 


real 


r.c. 


l.c. 


rest 


AO 


0.04 


0.04 


0.00 




-0.04 


-1.04 


-0.08 


1.08 


BO 


0.04 


0.03 




0.01 


-0.04 


-1.15 




1.11 


CO 


0.04 


0.03 




0.00 


-0.04 


-1.93 




1.89 


A 


0.20 


0.05 


0.15 




-0.20 


-1.06 


0.86 




B 


0.20 


0.03 


0.17 


0.00 


-0.20 


-1.02 


0.81 


0.01 


C 


0.20 


0.03 


0.17 




-0.20 


-1.18 


0.97 





Table 2: The ^-channel p -meson contributions to the polarizabilities a + /3 and a — j3 for various models 
described in the text. The column "real" lists the values obtained directly from the real parts of the model 
amplitudes, the columns "r.c." and "l.c." show the contributions of the dispersion integrals over the right and 
left cuts, respectively, the column "rest" gives contributions from Cauchy integrals around unphysical poles 
and at infinity. Within rounding errors, the three integral contributions add up to the entry "real". 

Unfortunately, also somewhat modified resonance models fail to give reasonable predictions based 
on dispersion relations. In Ref. [0], we have studied 6 such models: Model A with a constant width 
F(s) = To, model B with a pole at M — iF(s)/2, and model C derived from model B by dropping 
terms quadratic in T(s) as well as models AO, B0, and CO related to the previous models by in- 
troducing an energy-independent coupling, g(s) = g(M 2 ). As shown in Table ^J, all these models 
contain contributions from unphysical cuts, poles, or the "big circle" at infinity 

4. Conclusions 

The polarizabilities of the pion are elementary structure constants and therefore fundamental 
benchmarks for our understanding of QCD in the confinement region. These polarizabilities have 
been calculated in ChPT to the two-loop order with an estimated error of less than 20 % [g]. It 
is therefore disturbing that these predictions are at variance with the results of L.V. Fil'kov and 
V.L. Kashevarov by 2 standard deviations. The latter results are based on meson exchange models 
and performed in the framework of dispersion relations. As an example, this procedure leads to 
a?r+-/V = 13.60±2.15 [§, whereas ChPT predicts 5.7 ± 1.0 [g], all in units of 10~ 4 fm 3 . The 
discrepancy originates from huge contributions of intermediate meson states in the approach of 
Ref. [Q|. In ChPT, on the other hand, the vector mesons enter only at 0(p 6 ) through vector-meson 
saturation of low-energy constants. They are usually treated in the zero-width approximation and 
estimated to yield a much smaller effect, e.g., an ft) contribution of about 1 unit to the neutral 
pion polarizability [§, §j. Moreover, Fil'kov and Kashevarov predict large contributions of vector 
mesons to both the electric and the magnetic polarizabilities of the pion, at variance with the fact 
that the transition from the pion (J p = 0~) to the intermediate vector meson (J p = 1 ) is driven by 
magnetic dipole radiation leading to a (para)magnetic contribution only. 

The apparent discrepancy between the two approaches can be traced to the specific forms for the 
imaginary part of the Compton amplitudes [Qj, which serve as input to determine the polarizabili- 
ties at the Compton threshold (s = m 2 , t = 0) by dispersion integrals. In order to obtain amplitudes 
with good properties at high energies, Fil'kov and Kashevarov introduce energy-dependent cou- 
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pling constants with a square-root singularity (e.g., 1 / y/s or 1/ The resulting amplitudes fulfill 
the following conditions to set up dispersion relations: 

(I) The amplitudes are analytic on the physical Riemann sheet except for isolated points on the 
real axis. As an example, the ^-channel singularities of Ref. [^] are situated (i) at the threshold for 
two-pion states (s = Am 2 ), which leads to the physical cut, and (ii) at the origin of the Mandelstam 
plane (s = 0), which leads to an unphysical cut. 

(II) The amplitudes are square integrable along any line parallel to the real axis, albeit at the ex- 
pense of the unphysical cut starting at s = 0. 

Based on these conditions, Titchmarsh's theorem asserts that the real and the imaginary parts of 
the amplitudes are Hilbert transforms, that is, they are related by dispersion relations. 
In their actual calculations, Fil'kov and Kashevarov ignore the contribution from the unphysical cut 
by setting the imaginary part equal to zero "below the threshold of two-pion production". However 
one implements this statement in practice, a look at the contour plot in Fig. |3| shows that the proce- 
dure will inevitably introduce a "wall" in the complex plane, separating a region with finite values 
of the imaginary part from zero values. Because all the points on this wall become non-analytic 
ones, the function itself is no longer analytic and therefore condition (I) is no longer fulfilled (note: 
condition (I) allows only for isolated singularities). Furthermore, in the region of vanishing imag- 
inary parts, the only possible analytic function is a real constant everywhere. And unless this 
constant vanishes, one easily finds that also condition (II) breaks down. These arguments show 
the inconsistency of first introducing the factor 1 / \/s for convergence at large energies and later 
ignoring its consequence, the unphysical cut, at low energies. The results of this inconsistency 
are: (i) the vector meson effects are grossly overestimated and (ii) the magnetic dipole transition 
7 +71 1> p/co yields a large negative contribution to the electric polarizability. Similar problems 
show up for the exchange of other mesons. In particular, the l/y/t factor for a exchange in the t 
channel leads to a diverging amplitude M ++ (t) at t = 0, the point at which the polarizability is to 
be predicted. Even apart from dispersion relations, we would not recommend to fit the data in the 
measurable region by functions that approach infinity at or near the point to which one wants to 
extrapolate. In conclusion, the reported discrepancies between ChPT and dispersion theory result 
from applying the latter theory to non-analytic functions. 
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